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Abstract 



We show that it is possible to formulate Abelian Chern- Simons theory on a lattice 

00 

as a topological field theory. We discuss the relationship between gauge invariance of the 

O ' 

q ■ Chern-Simons lattice action and the topological interpretation of the canonical structure. 
(N 

We show that these theories are exactly solvable and have the same degrees of freedom as 
o ■ the analogous continuum theories. 

> 

X 
=3. 

In the continuum, Chern-Simons theory [1,2] is a topological field theory [3,4]. It 
is exactly solvable and is related to interesting topological structures such as the Witten 
invariants of three-manifolds and the Jones polynomial and other knot invariants for links 
embedded in three-manifolds. Its canonical formalism also has an interesting relationship 
with conformal field theory in two dimensions. As a model for physical phenomena its U(l) 
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(and sometimes U(N)) versions are proposed as the infrared limit of phenomenological 
theories of the fractionally quantized Hall effect [5] and other theories of anyons, including 
speculations about the mechanism for high T c superconductivity [6]. 

In recent years there have been efforts to write Chern-Simons theory on a lattice [7- 
11]. This gives an automatic ultraviolet regularization of the field theory and is important 
for condensed matter and statistical physics applications. It also makes well-defined the 
correspondence between matter-coupled Chern-Simons theory and a theory of anyons [11]. 
Generally, it is difficult to formulate lattice Chern-Simons theory in a natural way as one 
lacks geometrical principles which would replace the general covariance of a topological 
field theory in the continuum. 

Nevertheless, we shall show in this Letter that when lattice Chern-Simons theory is 
formulated so that it is both local and gauge invariant on the lattice (here the term "local" 
means that the canonical structure is such that a link commutes nontrivially only with 
links with whom it shares a common site), we obtain an exactly solvable model which 
shares many of the features of its continuum kin and can be regarded as a topological field 
theory on the lattice. 

To begin, we review some of the features of continuum U(l) Chern-Simons theory 
on the space M s = R 1 ® E s , where E 9 is a compact oriented two-dimensional Riemann 
surface with genus g. The action is 



The vector density j M (x) = J dr^i ^ r i( T )^ 3 ( x ~ r i( T )) corresponds to a collection of 
Wilson loops parametrized by rf (r) and <9 M j M = 0. As is usual in gauge theory, the action 
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(1) 



is linear in the temporal component of the gauge field and the field equation resulting from 
its variation enforces the constraint 

B(x) + T f(x) = (2) 

where B = dA (restricted to E 3 ) is the magnetic flux. It is useful to think of the remaining 
first order action, f dt J Es ^— -^A A A — A ■ j^j as already being cast in phase space, which 
is the set of gauge connection 1-forms on E s . From this we obtain the canonical commutator 
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Ai(x),Aj(y) = — e 0ij S 2 (x - y) (3) 

which has the property that, if we consider two distinct oriented curves, C and C on E s , 

= ?fv[C,C'] (4) 



A, I A 



*c jc 

where z/[C, C] is the number of right-handed minus the number of left-handed intersections 
of C and C (A right-handed intersection occurs when, if we move along the positive 
direction of C, C crosses from right to left.) Thus, the commutator gives a representation of 
the intersection form for Wilson loops. For closed curves C and C", v[C, C] is a topological 
invariant and if either C or C is contractible, v[C, C] must vanish. Therefore, there is a 
nontrivial commutator only for homologically nontrivial curves. If otj, /3j, j = 1, . . . ,g, are 
a canonical set of closed curves on E 9 generating its first homology group then v[cti, ctj] = 
v\fli, Pj] = 0, v[oti, (3j] = Sij. Upon imposition of the constraint (2) the integrals of A over 
the 1-cycles oti, fa are the only remaining degrees of freedom and, modulo the remaining 
symmetry under large gauge transformations, they form the reduced phase space [2]. 

It is the property (4), that the symplectic structure gives the intersection form for 
loops on E 3 , which we shall discover on the lattice, as a consequence of lattice gauge 
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invariance and locality. We work in continuum time and a square spatial lattice with 
spacing 1. We begin by fixing some notation. The forward and backward shift operators 
are Sif(x) = f(x + i) , S~ 1 f(x) = f(x — i) , respectively, and forward and backward 
difference operators are dif(x) = f(x + i) — f(x), di = Si — 1, dif(x) = f(x) — f(x — 
i) , di = 1 — S^ 1 = S~ 1 di. Summation by parts on a lattice takes the form (neglecting 
surface terms) J2 X f(x)dig(x) = —J2 x dif(x) g(x) by virtue of the lattice Leibniz rule 
di(fg) = fdig + difSig = fdig + Si(difg) (no sum on i). The components Ai(x) of the 
gauge field are real- valued functions on the links specified by the pair [x,i], A is a function 
on lattice sites, the magnetic field B(x) = diA 2 (x) — dzA^x) is a function on plaquettes 
where x labels the plaquette with corners x,x+l,x+l + 2,a; + 2, and the electric field 
F 0i = Ai — diA is a function on links. 

A gauge invariant, local, nondegenerate Chern-Simons term was found in ref. [11] : 

S= j dtJ2(^Mx,t)e lJ d l A J (x,t) - ^-Aiix^KijAjfat) + A^x,t)f(x,t)\ (5) 

where i,j = 1, 2, with a conserved current, dojo — diji = 0^ and 

K __\( s 2 -s^ - { - 1 + s - 1 + s 1 + s 2 - 1 s l ] 

13 2\-l + S^ +S 2 + S^S 2 S^-S! 
1 ( d 2 + d 2 -2 - 2di + 2d 2 + d 2 d x 



2 V 2 + 2d 2 - 2d! - d x d 2 -d x - d\ 

Analogous to (3), determines the symplectic structure on the phase space which is 
the space of functions Ai(x) from the links of the lattice to R 1 , 
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[A i (x),A j (y)] = -—K- j \x-y) (7) 



^ The divergence of a vector field is correctly transcribed onto the lattice through back- 
ward differencing. 



Also, the equation of motion for A gives the gauge constraint 
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e^d i A j (x) + —j°(x) = 



(8) 



This associates the magnetic flux in the plaquette with corners x,x+l,x+l + %x + 2 
with the charge at x. It was shown in [11] that (7) gives Wilson line-sums the lattice 
analog of the algebra (4), 



^2 dxiAi(x) 7 ^2 d yj A j(y) 

[x,i]eC [y,j]eC 



2ni 



(9) 



where the sums are over lattice curves (connected assemblies of oriented links) and dxi = 
±1 according to whether the link [x,i] is backward or forward directed. 

We may work backwards from (9) , and begin by considering it a primitive requirement 
that K~ x is anti-Hermitean, nondegenerate and gives the symplectic structure in (9). In 
particular, if C in (9) is a homologically trivial closed loop and C is an open curve, 
then v[C,C] is an integer, the number of times C links each endpoint of C", with sign 
determined by the orientations of C and C . This means that K~ x must obey 



di = dj6jkK ki ,di = —K^ ejkdk , 



(10) 



This turns out to be identical to the condition on K that one obtains by requiring that 
the action (5) is invariant under the gauge transform A^ — > Ai + diX, A — > A + x, or 



equivalently that -J^B(x) is the generator of static gauge transforms, 



k 

2tt 



= -idiXiv) 



(11) 



Thus, the local form of the gauge constraint in (8) as well as gauge invariance of the action 
are equivalent conditions to the topological invariance and integer- valuedness of v[C,C] 
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in (9) when either C or C is a closed contour which is the boundary of a set of plaquettes. 
Pictorially, gauge invariance, together with the locality of the symplectic structure K~ x 
guarantee that curves which touch but do not penetrate are counted as zero intersection: 



Conversely, these conditions, together with the requirement that K allows interac- 
tions of links with only those other links that share a common site, and that it have a local 
inverse, may be summarized by precisely the condition (10) (which we found by requir- 
ing gauge invariance), and makes K~ x the kernel that counts intersections of two lattice 
curves. The remaining structure of K~ x is the normalization, and is fixed by requiring 
that (9) should hold for the simple crossings 
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This, together with the Gauss Law" 1 " determines K 1 and K uniquely as that given in 
equation (6). Several results follow for partial intersection numbers: 



These would be ambiguous in the continuum and on the lattice they follow uniquely 
from our other requirements on K. 

Note that the Kernel in the action, Scs = / dt ^ f i(x)D fIU (x — y)A u {y) here can be 
parametrized as D^ u = T^\t\ au d CF and that the requirement of gauge invariance is precisely 
dT = d (where do = do = do), the analog of (10). When we require that T is a local matrix, 
with a local inverse, (the determinant is a monomial in shifts) this T _1 -matrix is the kernel 
which counts the intersections of surfaces with curves on a 3-dimensional lattice. This, in 
turn implies that (as we shall see in the following) the effective action for Chern-Simons 
theory coupled to currents is related to the linking numbers of trajectories. 

To solve the model (3), we first consider a lattice with trivial first homology group, 

t In fact, there are exactly four such if's, corresponding to the four possible ways of 
assigning a plaquette to one of its corners. Gauss' Law in these cases would appear as e.g. 
B + ^Sdo = 0. 

7 



such as the latticized open plane R 2 . We represent the commutator (5) by choosing the 
functional variables 

„, . 2ni d w . if-, dK~ 1 d „, . \ 

B « = -TaW) • A <*> = ^ ( rf • A - "^T B <*> J < 12 » 

Then, a wave-function which solves the constraint (8) is 

*phy 8 [A,^,t]=exp^i^A(x)j (x,t)j (13) 

The Hamiltonian is 

3 ^ V d-d d-d / 

where we have used the identity 5y = eikdk-j-^ejidi + di-^~dj. The Schrodinger equation 

^*phys[A,j M ,t] =H* phya [\J tl ,t] (15) 

is solved by 

(16) 

The trajectory of a particle is a piecewise linear lattice curve consisting of instan- 
taneous hoppings in spatial directions between lattice sites and temporal segments rep- 
resenting the particle at rest on a particular site. The second term in the phase of the 
wave-function in (16) is a topological invariant and can be interpreted as the linking num- 
ber of periodic lattice trajectories. To see this, first consider adding a closed spacelike 
curve to j M , the perimeter of a plaquette. This is described by the change in current 

Sj°(x,t) = , Sj 1 (x,t) = -d 2 S(x-a)S(t) , 5j 2 (x, t) = di5(x - a)5(t) 
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which has the property that d x Sj(x,t) = d ■ dS(x — a)S(t), implying that the phase of 
the wavefunction changes by ^-9(t)j°(a, 0), which is times the linking number of the 
lattice curves with the plaquette. If we add a time-like plaquette, for example 

5j°(x,t) =e(t-t l )6(t 2 -t)d 1 5(x-a) , Sj\x,t) = (5(t - t x ) - S(t - t 2 ))S(x - a) 

The phase changes by 

— jf dtj 2 (a,t) 

which is also 2tt jk times an integer, the number of times the lattice curves link the time- 
like plaquette. To get the latter result we have neglected the 'self-linking' of the space- 
like plaquette with itself. Such self-linking numbers are not well-defined here but require 
further regularization. (Here, self-linking numbers involve ill-defined products of the form 
0(t)6(t).) 

In general, a lattice current due to a charged particle which hops from lattice site x v -\ 
to x p along links C p at time t p has the form 

f(x, t) = J2 6 ( x - z P )0(Vn - t)e(t - t p ) 

V 

j\x,t) = Y J E S(x-y p )dy ip S(t-t p ) (17) 

P [VpAp] =links in C p 

A general conserved current may be written as a sum of currents of the form (17) = 
J2 a j%' quadratic self-interaction terms for currents in (16) are of necessity ambigu- 
ous, because they are in fact the self-linking number of the trajectories - the latter is a 
topological invariant of framed links, and the ambiguity can only be lifted by introducing 
a framing. 
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The cross terms of (16) between two trajectories a, (3 are also ambiguous unless the 
instantaneous hoppings in spatial directions of the two lattice curves occur at distinct 
times. When this condition is met, the cross term is well defined. The a, f3 cross term may 
be combined with the /3, a cross term, and calculated in terms of a lattice angle function. 
Lattice angle functions have been discussed in earlier works [8,11], and are usually defined 
as a contour sum as 9c{x,y) = 2nY^d£i^-*8{i ~ v) an< ^ are multivalued functions of 
position in that they depend on windings of the contour C around the point y. 

d x d6c{x,y) = 2ttS(x - y) (18a) 
Oc(x,y) - 6c>(x,y) = 27rJ2^(CC'-\z) (186) 

z 

where w(CC" _1 , z) is the winding number of the closed curve CC'~ X around the point z. 
Unlike the continuum angle function, which satisifies the further identity 

6 c (x, y) - 9 C > (y, x) = n + 2nis[C, C'\ (19a) 

this angle-function has the property that 

6 c (x, y) - 9c (y, x) = n + 2nu[C, C] + e(x - y) (1%) 

The last term in (16) is precisely the defect function £(# — y) for the lattice angle function 
which the first term in (16) calculates. Thus the result for the cross terms a, (3 and /?, a 
together is proportional to the total change in angle between the curves accumulated during 
the time evolution, 

1 f d 

- I dt -^{x{t) ~ x'(t)) , x ea x G f3 

(Note that the result is independent of C.) where 9 is an "improved" lattice angle function, 
in that it satisfies (19a) as well as (18ab). Thus for a periodic trajectory of iV particles 
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on the lattice the phase of the wave-function is the linking number of trajectories and the 
wave-function therefore carries a 1-dimensional unitary representation of the iV th -order 
braid group of the plane where braiding is constrained to follow links of the lattice. This 
implies that particles coupled to the lattice Chern-Simons theory are anyons with statistics 
parameter 1/k. Furthermore, aside from this phase the theory is trivial. The Hilbert space 
contains only one state. 

A more complicated situation arises when the lattice has a nontrivial first homology 
group. Here, for simplicity we shall set j M = and consider the example of a toroidal 
lattice where the gauge fields have periodic boundary conditions 

At(xi + iVi,x 2 ) = A i (x 1 ,x 2 ) , A i (x 1 ,x 2 + N 2 ) = A i (x 1 ,x 2 ) (19) 

Here, the gauge group is U(l) and we require that the gauge transformation obeys the 
boundary condition 

X(xi + N u x 2 ) = x(xi,X2) + 27rmi , xi x U x 2 + N 2 ) = x(xi, x 2 ) + 2Tvm 2 (20) 

where are integers. We can choose the canonical generators of the first homology and 
the gauge invariant canonical variables as 

q = J2Mn,0) , p = -|^$> 2 (0,n) (21) 

n n 

which obey the algebra 

[q,p]=i (22) 

The large gauge transforms q — > q + 2irm\, p — > p — km 2 are generated by mi and m 2 
operations of the unitary operators g\ = e 2lTip g 2 = e lkq , respectively. These operators 
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obey the algebra 

9192 = 929ie 27rtk (23) 

The wave-function should carry a unitary represention of this algebra. Here we are assum- 
ing that k = k\/k 2 is a rational number. It is straightforward to construct a representation 
of (23). Assume that we find an eigenvector ipg of g\ such that 

g^e = • (24) 

Operating g 2 £ times gives another eigenvector of g\ with eigenvalue e l ( e + 27rki )^ 

9i 9^0 = 9 e 29ie 2 * ik % = e e+ ^ iM 9^0 (25) 

When £ = k 2 we obtain the original eigenvalue. Thus we conclude that the algebra (23) 
may be represented by k 2 x k 2 unitary matrices and that g 2 2 ipg^ = e^tfjg^. Thus, the 
representation is specified by two angles 9 and £. In the Schrodinger polarization the 
wavefunction is a function of q and (24) implies 

M(q) = J2 ein+6/27Z)iq ^dn) (26) 
nez 

and 

<&M«) = E e^ +0 ^ +ki ^^(n) (27) 
nez 

When £ = k 2 we obtain ip0^(n) = e^ipe^(n + ki). Thus, for each £ there are k\ independent 
coefficients in the expansion (26). The Hilbert space is /c^-dimensional. This construction 
can be generalized to lattices with more complicated homology. For a lattice with genus 
<7, the dimension of the Hilbert space turns out to be (kik 2 ) 9 . This dimension is the same 
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as that of the Hilbert space in the continuum theory on the space R 1 x E s where E 3 is a 
Riemann surface of genus g. 

In conclusion, we have shown that the present version of lattice Chern-Simons theory, 
where the form of the action is deduced from the requirements of locality of the action 
and canonical structure, as well as from gauge invariance, recovers the topological features 
of continuum Chern-Simons theory. An interesting feature is the introduction of a new 
and practically unique lattice kernel which can be regarded as an intersection density and 
which is responsible for the topological nature of the theory. We showed that the lattice 
model is exactly solvable. The solution resembles the solution of the analogous continuum 
Chern-Simons theory. 

The wave-functions carry a representation of the Braid group. Note that, since par- 
ticles are constrained to occupy lattice sites, their motion is more restricted than in the 
continuum, and it is possible to form clusters of particles such that not every braiding 
operation is allowed. Thus, the braid group on the lattice must differ from the one in the 
continuum by certain constraints. It would be interesting to investigate this further. 

It is intriguing to us that the Hilbert space of lattice Chern-Simons theory on the torus 
is identical to that in the continuum. In previous work on the continuum theory [2] the 
Hilbert space was related to quantization of the moduli space of the torus, i.e. the space of 
metrics on the torus modulo diffeomorphisms. On the lattice there are no diffeomorphisms 
or metrics and no concept of moduli space. Yet we find a quantum theory on the lattice 
virtually identical to that in the continuum. 
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